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ABSTRACT 

We investigate the narrow-line region (NLR) of two radio-quiet QSOs, PG1012H-008 and PG1307H-085, us- 
ing high signal-to-noise spatially resolved long-slit spectra obtained with FORSl at the Very Large Telescope. 
Although the emission is dominated by the point-spread function of the nuclear source, we are able to detect 
extended NLR emission out to several kpc scales in both QSOs by subtracting the scaled central spectrum from 
outer spectra. In contrast to the nuclear spectrum, which shows a prominent blue wing and a broad line profile 
of the [O III] line, the extended emission reveals no clear signs of large scale outflows. Exploiting the wide 
wavelength range, we determine the radial change of the gas properties in the NLR, i.e., gas temperature, den- 
sity, and ionization parameter, and compare them with those of Seyfert galaxies and type-II QSOs. The QSOs 
have higher nuclear temperature and lower electron density than Seyferts, but show no significant difference 
compared to type-II QSOs, while the ionization parameter decreases with radial distance, similar to the case of 
Seyfert galaxies. For PG1012-I-008, we determine the stellar velocity dispersion of the host galaxy. Combined 
with the black hole mass, we find that the luminous radio-quiet QSO follows the local Mbh - relation of 
active galactic nuclei. 
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1. INTRODUCTION 

Since the discovery of correlations between the mass 
of black holes (BHs) and the prope rties of their host- 
galax y bulges (Magorrian et al. 1998; iFerrarese & Merrittl 
120001; iGebhardt et al., ,2000.) . the interplay between active 
galactic nuclei (AGNs) — thought to represent a stage of 
galaxy evolution in which BH is actively growing through 
accretion — and their host galaxies has received much at- 
tention. In particular, AGN feedback seems to provide a 
promising way to drive these relations by quenching both 



BH growth (e.g., 


Di Matteo et alj|2005l; iHopkins et alj|2006l; 


ISironi & Socrates 


I20I0I). 



Recently, there has been growing observational support 
for such a s cenario for various AGNs, e.g., radio-loud 
quasars (e.g., iNesvadba et all 1200 8; Fu & Stockton 20090, 
broad-absorptio n line quasars (e.g., .Crenshaw et al.. ,20031 
Moe et al1l2()09D. local ultra-luminous inf r ared galaxies (e.g. , 



and narrow-line region (NLR), which are photoionized by the 
central engine. While BLRs are confined to sub-pc scales 
around an accretion disk, producing kinematically broadened 
emission lines with typical velocities of lO-' - lO'* kms~', 
NLRs, characterized by narrow-lines with typical velocities 
of 10^- 10^ kms"', can span over kpc scales, which are com- 
parable to the s ize of the bulge or even the entire galaxy 
(iBoroson & Ok5 1984; Stockton & MacKentvl [19871) . Thus, 
as a direct interface between the AGN and its host galaxy, the 
NLR can provide important clues on the impact of BHs on 
their host galaxies and vice versa. 

The NLR of Seyfert galaxies has been studied ex tensively 
over t he last decades (e-g - . iM ulchaey et al. 199 61 iHo et al 



1997 



20021 



Fischer et al.' '2010; 'Feru glio et all l20Tot ISfiirm et al.l I201U 
Rupke & Veilleux 2 011), and high redshif t AGNs (e.g. 
Tremonti et all 120071: [Alexander etaTI 120 lOt iHainline et d] 
201 Ih . These observations are, however, often restricted to 
a handful of sources. Moreover, it is difficult to determine 
the kinetic energy involved and to prove that the observed 
outflows are indeed galaxy-scale radiatively-driven AGN out- 
flows (compared to outflows driven by e.g., star formation or 
jets). 

One approach in the quest for observational signatures of 
AGN feedback is to focus on the emission-line regions in 
the vicinity of the BH: the so-called broad-line region (BLR) 

' Author to whom correspondence should be addressed. 



Schmittetal. 2003; Fal cke et al.nl 998c iBennert et al 
Stoklasova et al. 2009). While most earlier studies 
have focused on the extent and morphology of the NLR, 
the interplay between radio emission and N LR, and implica- 
tions for the unified model of AGNs (see lAntonuccil [l993l 
for review), the goal of recent studies is to search for im- 

J Tints of AGN feedbac k, e.g., in the form of o utflowing gas 
Schlesingeretal ."2009; Mulle r-Sanchez et al ."201 1*) . 

In particular, the forbidden [O III] A5007 (hereafter [O III]) 
emission line is well known to show blue w ings, which are 
gener ally interpreted as a sign of outflow (see ICrenshaw et all 
120101 and references therein). To study the kinematics of 
the ionized gas, spectroscopy of the extended NLR is essen- 
tial. However, such a study is especially challenging for the 
most promising candidates to exhibit AGN feedback, high- 
luminosity QSOs, due to the presence of their bright nuclei. 

Recent studies have thus focused on type-II (obscured) 
QSOs that have been discovered in large numbers from 
the Sloan Digital Sky Survey (SDSS; .Zakamska et al. 20031: 
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Figure 1. Central de-redshifted spectra of PG1012+008 (top) and 
PG1307+085 (bottom). 



iReves et alj 120081) . For example, iGreene et all (120111) stud- 
ied 15 luminous type-11 QSOs at low-redshift (z < 0.5) with 
spatially resolved spectroscopy. Comparing the extent of the 
continuum and [O III] emission, they argue that the AGN is 
ionizing the interstellar medium in the entire host galaxy. The 
large velocity dispersion of this ionized gas out to kpc scales 
suggests that the gas is disturbed on galactic scale. However, 
the velocity dispersion is below 500 kms"' for most objects 
in their sample, and the overall escaping fraction is less than 
25% with a median of 2%, e v en in the most extreme cases. 
Similarly, Villar-Ma rtm et al.l (1201 ll) report blue asymmetry 
of the [O III] emission line in 1 1 objects out of a sample of 13 
type-11 quasars at 0.3 < z < 0.5. They interpret the asymmetry 
as a signature of outflows within a few kpc from the central 
engine. Furthermore, they find an anticorrelation between the 
degree by which the kinematics of the outflow is perturbed 
and its contribution to the total [O III] flux, suggesting that 
only a small fraction of the total mass of the ionized gas is 
involved in the outflow. 

With this line of evidence for type-11 quasars, we here study 
the NLR of type-1 QSOs, PG 1 1 2H-008 and PG 1 307H-085 with 
spatially resolved spectroscopy. Benefiting from the wide 
spectral range and excellent seeing condition, we investigate 
the kinematics and physical conditions of the NLR as a func- 
tion of radial distance. Throughout the paper, we assume a 
Hubble constant of //q = 70 km s~' Mpc"', fl^ = 0.7, and I^m 
= 0.3. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Sample Selection 

To study the NLR of luminous quasars, we started 
with the seven low-red shift Palomar-Cireen (PG) quasars 
dSchmidt & GreenI Il983h . for which [O III] narrow-band 
imaging was obtained with the Hubble Space Telescope (HST) 



by iBennert et al.l (120021) . Out of these seven objects, we se- 
lected those with the most extended NLR that were visible 
from the Very Large Telescope (VLT), and at a redshift range 
that allowed us to cover the broadest range of emission lines 
in optical spectra. As a result, we selected two such objects, 
PG1012H-008 (z = 0.187) and PG1307H-085 (z = 0.155), with 
an NLR size of 1 "59 (- 5 kpc) an d 1 "39 (~ 3.7 kpc), r espec- 
tively, based on the HST images (' Bennert et al.ll20()2h . The 
morpholog y of the [O III] e mission is symmetric and compact 
in general (Leip ski et al.ll2006i) . The surface brightness at the 



1.4 X 10"'^ erg s' 



"'cm ^arcsec ^ 



(3-cr 



edge of the NLRs is ' 
detection). 

These two QSOs wer e also studied with the Wide-Field 
Camera on-board HST dBahcall et alj 119971). an d with the 
Very-Large Array in the radio (iLeipski et al.ll2006b . Here, we 
briefly summarize the properties of the host galaxies and the 
radio emission of the sample based on these previous studies. 
The optical image of PG1012H-008 reveals a complex mor- 
phology: It is undergoing a major merger with a companion 
galaxy at 3"3 10 kpc) to the east, exhibiting clear tidal-tail 
structure. Another smaller, compact galaxy at 6. "8 (^21 kpc) 
north of the QSO may also be interacting with the two galax- 
ies. In the case of PG1307H-085, the host galaxy is an isolated 
small elliptical galaxy. 

Both QSOs are radio-quiet with total fluxes less than 
0.3 mJy at 8 .4 GHz for PG10 12-H008 and 4.8 GHz for 
PG1307H-085 (Leipski et ajj|2006|). The isophotes of [O III] 
and the radio continuum show similar structure in both ob- 
jects on a scale of ^ 1", implying a radio-NLR interaction. 
Particularly for PG1012H-008, a "bending" in the isophotes is 
apparent in the [O III] and radio images in the SE direction, 
although it is unclear what role is played by the merger 

The radial change of the [O III] line profile of both PG 
QSOs and its connection to the radio emission was investi- 
gated by Leipski & Bennert (2006) based on the spectra ob- 
tained at the New Technology Telescope at the European 
Southern Observatory (ESO). Investigating the substructures 
of the [O III] line profile, they note that jet-gas interaction 
plays a role in the NLRs of radio-quiet quasars. Although 
their method of analysis, i.e., spatially-resolved spectroscopy, 
is similar to ours, we present a more complete view on 
the NLRs using superior data with better spatial resolutions, 
wider wavelength ranges, and far higher signal-to-noise ra- 
tios (S/Ns), though we miss some kinematic subcomponents 
identified due to relatively lower spectral resolution. 

2.2. Observations 

To investigate the properties of the NLR based on spatially 
resolved specti-a, we obsei-yed PG1012H-008 and PG1307H-085 
with the visual and near-UV FOcal Reducer and low disper- 
sion Spectrograph (FORSl) at the VLT on 2005 April 1 and 
3. To achieve high spatial resolution, the observations were 
carried out under good seeing conditions (< 0"7). A slit 
width of 0."7 was chosen to match the seeing. Based on the 
[O III] narrow-band images, we aligned the position angle of 
the slit with the maximum radial extent of the [O III] emis- 
sion. The spatial resolution is "2 pixel"'. We used grism 
300V with an order sorting filter GG375, covering a spectral 
range of 3300-7180 A with 2.6 A pixel"'. The instrumental 
broadening (uinst), determined from arc lamps and sky lines, 
is 216 kms"' at 5000 A. For flux calibrations, we observed 
a spectrophotometric standard star G60-54. Table [T] summa- 
rizes the observations. 
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Table 1 

Observation 



Name z Exposure P.A. Seeing S/N Scale Dl i[o//;j Lb Re 

(s) (deg) (") (pixer') (kpc arcsec"' ) (Mpc) (x lO'*^ erg s"') (x lO'*^ erg s"' ) (") 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) 

PG1012+008 0.187 7400 142 < 0.7 472-26 3.1 909.4 4.94 1.13 3.4 

PG1307+085 0.155 6800 112 < 0.6 507-5 2.7 739.2 4.53 1.41 1.3 

Note. — Column 1. object name. Column 2. redshift determined from [O II] A3727. Column 3. total exposure time. Column 4. position 
angle of the slit. Column 5. seeing condition. Column 6. signal-to-noise ratio at 5100 A from the center to the outermost aperture. Column 
7. angular diameter distance. Column 8. luminosity distance. Column 9. [O III] luminosity measured from Be nnert et al. (2002). Column 
10. B-band luminosity from Elvis et al. (1994). Column 11. effective radius derived from two-dimensional de Vaucouleurs fits to the surface 
brightness profile given in Bahcall et al. ( 1997). 



2.3. Data Reduction 

Standard data reduction, i.e., bias subtraction, flat-fielding, 
wavelength calibration, flux calibration, and atmospheric ab- 
sorption correction was carried out using a serie s of IRAP 
script s dev eloped for long-slit spectroscopy (e.g., IWoo et al] 

f35 , 2006). Cosmic rays were removed using L. A. Cosmic 
nDokkum 2001). We performed wavelength calibra- 
tion using arc lines and obtained a spectral resolution of 
2.6 A pixel"' with standard deviation of ^ 0.25 A. 

We extract seven spectra along the slit (spatial direction), 
covering ±3 "6. Positive direction indicates SE and NW for 
PG1012H-008 and PG1307H-085, respectively. The aperture 
size was increased from 0."6 at the center to 1 "8 at the outer 
radii by a step of 0."4, in order to boost the S/N. Figure [T] 
presents the spectra extracted from the central aperture for 
both targets. Note that the region near Ha suffers from the 
atmospheric absorption (A band). We attempted a correction 
using the standard star spectra. However, we were unable to 
measure the [N II] and narrow Ha emission lines, due to the 
dominance of the broad Ha line. 

3. RESULTS 

3.1. Extended Narrow-line Region 

We analyze spatially resolved spectra to investigate the 
spatial extension of the NLR. First, to compare the relative 
strength of the emission between the central and extended re- 
gions, we determine the spatial profiles of [O III], continuum, 
and standard star, respectively. In the case of the continuum, 
which is dominated by the AGN and thus expected to be spa- 
tially broadened by the seeing effect, we collapse the spec- 
tra in the continuum wavelength range 5080-5120 A to de- 
termine the spatial profile. For [O III], we sum the flux in a 
spectral range of ^ 30A centered on the peak of [O III], and 
subtract the continuum level previously determined. We also 
determine the point-spread function (PSF) based on the spec- 
tral image of the standard star. Figure 2 presents the spatial 
profiles of both targets, which are normalized by their maxi- 
mum flux value. In the case of PG1012H-008, the continuum 
spatial profile is consistent with the PSF in the inner region 
at -1" < r < 2", while at the outer region, the contribution 
of the host galaxy star light becomes dominant as indicated 
by a broad extension in the spatial profile. Although the spa- 
tial profile of [O III] is dominated by the PSF, the left side 
of the profile is slightly extended over the PSF or the contin- 
uum spatial profile. In the case of PG1307H-085, we detect an 
asymmetric profile of [O III] with a distinct bump at -3", rep- 
resenting an extended NLR. Thus, we conclude that the [O III] 
NLR in both objects is marginally detected. 



F — I 1 1 1 1 — 3 

\ continuum \ 




distance from the center (") 

Figure 2. Normalized spatial profile of continuum, [O III] emission line, and 
the standard star, respectively for PG1012+008 (top) and PG1307+085 (bot- 
tom). For the continuum, the shaded region indicates 3(t errors. 



To further test whether the [O III] emission in the off-center 
region is consistent with being the PSF wings of the [O III] 
emission from the center, we use the central spectrum of each 
object as a template to subtract the seeing-scattered broad and 
narrow lines from the outer spectra. For generating a tem- 
plate of the central spectrum, we fit the continuum and the 
broad component of H/3 (see the top panel of Figure 3), then 
combine the best fit models with the profiles of the narrow 
H/3 and [O III] . To account for the contribution from the host 
galaxy star light, we also add a linear background level. 

Figure |3] presents the continuum-subtracted central spec- 
trum (top panel) and residual spectra of outer region after sub- 
tracting a scaled central spectrum (middle and bottom panels). 
The strong broad H/3 component observed at the outer spec- 
tra is well subtracted along with the contribution of the PSF 
wings of narrow emission lines from the center, leaving a flat 
continuum and residual narrow emission lines. These resid- 
ual [O III] and H/3 lines represent the extended NLR in both 
objects. As expected from the investigation of spatial pro- 
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Figure 3. Continuum-subtracted central spectrum and residuals after subtracting a scaled central spectrum from the outer spectra, respectively, for PG 10 12+008 
(left) and PG1307+085 (right). Blue and red lines correspond to the minus and plus side of the center along the slit, respectively, and the number on the upper 
left corresponds to the distance from the center in units of arcseconds. The best-fit models using several subcomponents are represented by dashed lines (see the 
text for details), the residuals are presented with an arbitrary offset. The [O III] emission line is enlarged with the wavelength axis converted to velocity in the 
light panels. The reference redshift is determined from [O II] A3727 line in the central spectrum. 



files (Figure |2]i, both objects show very different spectra as 
the distance from the center increases. For PG1012-I-008, the 
extended [O III] emission is detected out to r ^ 2" with simi- 
lar flux on either side, while PG 1307-1-085 shows the extended 
emission only on one side. The host galaxy starlight, e.g., H/3 
absorption, starts to show at r > 2" for PG1012H-008, while 
it is not visible in PG1307-I-085, presumably due to a higher 
contrast of the nuclear AGN luminosity to the host galaxy 
star light. Thus, we detect the extended NLR emission out 
to r > 2.0ib 0."7 (6 lb 2 kpc), wh ich is larger than the sizes 
measured bv lBennert et al.l (120021) since the depth of our ob- 
servation is roughly 20 times larger. 

The central spectrum of both targets shows an asymmetric 
[O III] line profile with prominent blue wings, which can be 
interpreted as outflows. To properly measure the velocity of 
[O III], we model the H/3-[0 III] region of the central spec- 
trum as follows: First, we use a power law -i- IZwI template 
from Boroson & GreenI (11992b to fit strong Fe II multiplets in 
the central spectrum of PG1012H-008. For PG1307H-085, the 
Fe II emission is an order of magnitude weaker (compared to 
the continuum) and considered negligible. Then, we model 
the [O III] line profile with two Gaussians along with broad 
and narrow H/3. We use Gauss-Hermite polynomials of order 
six for the broad component of H/?, and the ratio of the two 
[O III] lines is fixed to the theoretical value (1/3). The best-fit 
models of the continuum-subtracted central spectra are shown 
in the top panels of Figure |3] 

The "blue" component of [O III], presumably represent- 
ing the outflowing gas, is blueshifted from the reference red- 
shift by 1380 kms-i and 173 kms"' for PG1012H-008 and 
PG1307H-085, respectively. The full width at half maxi- 



mum (FWHM) velocity of the total line profile measured 
from the central spectra is 1050 kms"' for PG1012H-008 
and 540 kms"' for PG1307-I-085, respectively. In contrast, 
as shown in Figure 3, the extended [O III] emission is nar- 
rower than the [O III] in the central spectrum for both ob- 
jects, and most of the line flux is coming from low veloc- 
ity gas (V < 500 kms"'). This indicates that the outflowing 
component with V > 500 km s"' is unresolved at the spatial 
resolution of our current observation and confined to within 
0"7 scale (~2 kpc). Thus, the emission-line profile of the 
extended gas indicates that the extended gas is kinematically 
quiescent in both QSOs. We conclude that while we detect 
extended emission out to r « 2" (6 kpc), there are no signs of 
a large-scale(~ 10 kpc) outflow. 

We briefly compare our results with those by 
iLeipski & Bennerj (|2006). The red asymmetry towards 
the south of the nucleus in PG1012H-008 detected by 
Leipski & Bennerj (|2iD06) may correspond to the redshifted 
residual [O III] line at 0.'^ 8. We m iss some of the complex 
velocity structure seen by ILeipski & Bennert (2006) due to 
lower spectral resolution. As noted by iLeipski & BennertI 
(120061) . PG 1307-1-085 shows a dramatic change in line profile 
especially to the SE ("minus" in this study) side. The 
residual narrow line to the SE jaf PG13 07-i-C) 8 is un resolved 
at our spectral resolution. ILeipski & iSennertl (l2006i) report a 
velocity dispersion of ^ 69 km s"' for this component, which 
is consistent with our results. Whi le our results are in overall 
agreement with those presented in 'Leips ki & BennertI (12006ft . 
given the differences in spectral resolution and S/N, we find 
that the broad [O III] Une profile of PG1307H-085 towards the 
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Figure 4. Ionization parameter f/ as a function of radial distance for 
PG1012+008 (top) and PG1307+085 (bottom). 
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Figu res. PG10 12+008 on the Mbh -c* relation. Black data and fit are taken 
from lWoo et alj HOl O). Our target is marked as a red circle with the arrow 
indicating that the stellar velocity dispersion measured at ±2" is a lower limit 
to the true luminosity-weighted bulge velocity dispersion. 



south is consistent with being a PSF wing rather than true 
physical extension of outflowing gas. 

The [O III]/H/3 l ine ratio is a diagn ostic of the ion- 
ization mechanism ([Baldwin et alj Il981h . In the case of 
PG1012-I-008, we cannot rule out that the ionization of the 
extended gas is at least partially due to star formation. For 
PG1307H-085, however, the high (> 5) [O III]/H/3 ratio at 
r < -0."8 from the center suggests that the QSO is the main 
source of ionization. 



3.2. Physical Conditions 

In this section, we measure temperature, density, and ion- 
ization parameter of our targets to compare them with other 
AGNs. 

3.2.1. Temperature and Density 

We determined the temperature of the NLR ga s using the 
rOlII] (A4959 + A5007)/A4363 ratio as oudined in .Osterbrockl 
([T98S 



,,A4959 + A5007^ 
A4363 ' 



7.73exp(3.29 x lOVr) 
1+4.5 X 10-4(n,/ri/2) ' 



(1) 



where, n^ is electron density and T is gas temperature. With 
typical values of = lO-' - 10"* cm"^ and T = lO'^K, the de- 
pendence on «(, is negligible. 

For PG1307-I-085, we were able to measure temperature 
and density using [O III] and [S II] line ratios from the cen- 
ter towards the south. The derived temperature is 25500 ± 
1480K at the center, and decreases with radial distance to 
13,4 00 ± 11.700K at 3f'6 (^ ^^9.7 kpc). Using the same meth- 
ods, iBennert et al.l (|2006a"b') estimated the average tempera- 
ture of four Seyfert 1 and four Seyfert 2 galaxies in the nu- 
clear region as 33,590 ± 7070K a nd 14,470 ± 440K, re- 
spectively, while lGreene et all (1201 lb estimated a temperature 
range of T = 11,000 — 23, 000 K for 15 type-I I QSOs. Note 
that the measurements of the Seyfert galaxies in'B ennert et alJ 
(l2()06a b) included reddening c orrection determ ined from the 
Balmer decrement, while Gre ene et al.l (1201 lb and this work 
did not correct for reddening. Note also that the physical scale 
in each measurements differs: for the Seyfert galaxies it cor- 
responds to a few pc in the center, while for the type-II QSOs 
in Greene et al. (201 1), it corresponds to 3" fiber aperture of 
SDSS spectra, covering 11 kpc at the median redshift. Al- 
though a direct comparison is not straightforward, the type -I 
QSO PG 1307-1-085 reaches lower temperature at the center 
compared to Seyfert 1 galaxies and similar to type-II QSOs. 

The density was measured in a standard way from the [S II] 
A67 16/6731 ratio using IRAF temden task, correcting for 
the effect of temperature using the derived values. The de- 
rived density ranges between 584 and 169 cm"^, decreas- 
ing with radial distance (0"-2."0). These values are simi- 
lar to the mean density = 3 35 cm"^ of type-II QSOs es- 
timated by ' Greene et al.l (1201 II) . and slightly lower than the 
mean nuclear density of Seyfert galaxies: 1 070 ±180 cm~^ 
for Seyfert 1 galaxies (IBennert et aLll2006al) and 1 100 ±315 
cm"^ for Seyfert 2 galaxies (IBennert et al.ll2006bl) . However, 
one should note that the nuclear density varies in a wide range 
(~ 300-2500 crn"^) a mong Seyfert 1 and 2 galaxies studied 
by IBennert et al.l (12006 a.b ). and it is not yet clear if the dif- 
ference is significant. The higher and denser NLR of type-I 
QSOs can be a natural result of a dusty torus blocking the 
central region in type-II QSOs. 

3.2.2. Ionization Parameter 

The ionization parameter is defined by the ratio of ioniz- 
ing density to hydrogen density as t/ = Q{H) /AircP'ne, where 
Q(H) is the number of the ionizing photons emitted by the 
AGN per second and n^ is electron density. By measuring the 

^ http://stsdas.stsci.edu/nebular/temden.html 
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Table 2 

Black Hole Masses 



Name iboi/iEdd" logisioo logLsioo'' cth/S logMBn/M© a*" 

(ergs"') (ergs"') (kms"') (kms"') 

PG1012+008 0.16 44.3 45.0 1729 8.42 196±30 

PG1307+085 0.14 44.2 45.0 1885 8.50 

" The bolometric correction factor for L5100 is taken from lNetze3 J2009D . 
From Vestergaard & Pet erson (2006J . 
Measured at 2" from the center. 



emission-line flux ratio [O II] A3727/[0 III] A5 007 the ion- 
ization parameter can be empirically derived (iPenston et alj 
[1990). We estimate the ionization parameter as a function 
of radial distance for both targets assuming = 100 and 
lie = 1000, respectiv ely, based on the calculations provided by 
iKomossa & Schulzl ( fl997L ; Figure g]). 

The ionization parameter peaks at the center, and decreases 
with radius for both of our targets. Given that the density 
decreases radially as well, the actual decrease of the ionization 
parameter is slower than shown in either one of the curves 
(solid versus dashed); more accurately, it would start in the 
center on the dashed line and then move toward values on the 
solid line with increasing distance from the center 

Note that we did not apply any reddening correction to the 
line fluxes. If the radial decrease of the ionization parameter 
were entirely due to change in reddening, a simple calculation 
using the standard extinction curve shows that AE(B-V) ^ 
2 would be required between the center and the outermost 
apertures. Considering that the aver age global reddening of 
Seyfert 1 galaxies is E(B-V) 0.4 dBennert et al.ll2006ah . it 
is unlikely that the radial trend is purely caused by reddening, 
although there might be a small contribution from reddening . 

Qualitatively, the decrease of the ionization parameter with 
radius is consistent with the radial trend found in Seyfert 
galaxies dBennert et al.ll2006allbh . The average nuclear ioniza- 
tion parameter of the two targets is L^„^=iooocm-3 = 1-79 x 10"^, 
which is a pproximately five tim es larger than Seyfert 1 and 2 
galaxies dBennert et al.ll2006allbh . 

3.3. Black Hole Masses and Stellar-velocity Dispersion 

We estimate BH masses from the velo city and size of the 
BLR based on virial assumptions (e.g.. IWoo & Urryl 120021: 
iPark et ani2012l) . In practice, we use the single-epoch virial 
mass estimator using the AGN continuum luminosity at 5100 
A (L5100) and the line dispersion of the broad H/3 line au/^- 
From the nuclear spectra extracted with an aperture of 2" ra- 
dius, including > 99% of the light, we measure L5100 and aiip, 
after subt racting Fe II contamination. Compared to the val- 
ues from Vestergaard & Peterson ('2006'), our L5100 is 0.77 dex 
smaller probably due to slit losses. Therefore, we use Lsioo 
from Vesterga ard & PetersonU i2006). 

Assuming a virial factor of log f = 0.72 dWoo et al.ll2010l) . 
we calculate Mbh- Recently, iPark et akl (120 12h reported the 
systematic difference of the H/3 line profile between the 
single-epoch and the r ms spectra . Corr ecting for this effect 
using Equation (13) in P ark et al.l (120121) . the mass estimates 
of Mbh are decreased by ~ 0.18 dex ^able|2]i. Th e estimates 
are co nsistent with those given in Vestergaard & PetersonI 
d2006l) within the errors of single-epoch BH mass (~ 
0.46 dex). 

For PG1012-I-008, we can also measure stellar-velocity dis- 
persion since we clearly detect the host galaxy (at 2" from 



the center). Due to the combined effect of the contamination 
from the AGN continuum and emission lines, and the pres- 
ence of 4000 A break, we consider the measurement based 
on the Ca H-i-K absorption lines as uncertain. In the case of 
the G-band region, broad H7 affects the fitting. Thus, stellar- 
velocity dispersion is determined using the spectral range 
5030-5300, including several strong stellar lines such as the 
Mg lb triplet (5 172 A) and Fe (5270 A) ('e.g.. lWoo etal.l2006l 
"200^: ' Greene & H0II2OO6I) . 

Both measurements now enable us to place the 
PG1012H-008 on the Mbh -cr, relation in Figure |5] where the 
24 reverberation-mapp ed AG Ns and the best-fit Mbh - cr* re- 
lation from Woo et al. (120101) are overplotted for comparison. 
Note that the stellar velocity dispersion of PG1012-I-008 is 
measured at 2" from the center, and can thus be considered a 
lower limit for the luminosity-weighted velocity dispersion 
of the bulge since the v elocity dispersion generally increases 
toward the center (e.g.. lKang e t al. 2013). Even when taking 
into account an expected shift to the right, the object falls 
onto the local Mbh - cr, relation (Figure |5), indicating that at 
least this low-z high-luminosity QSO is following the same 
Mbh-ct* relation. 

4. DISCUSSION AND SUMMARY 

We investigated the kinematics and physical condition of 
the NLR in two type-I QSOs, PG1012H-008 and PG1307H-085, 
using spatially resolved spectroscopy. The emission from the 
NLRs of these two targets is dominated by a central point 
source with slight deviation from the expected PSFs. How- 
ever, modeling the H/3-[0 III] region of the outer spectra 
clearly reveals the spatially extended weak NLR out to r ?s 2" 
(6 kpc) from the center, spanning a large r extent than previ- 
ously r eported from HST [O III] images bv lLeipski & BennertI 
d2006l) due to the increase in depth of our observations. The 
[O III] line of the nuclear spectrum of both targets shows blue 
wings and broad widths with FWHM velocities exceeding 
500 kms"'. However, the same feature is not present in the 
outer spectra, indicating that the outflow is confined within 
the PSF size (~2 kpc) and does not extend to galactic scales. 

Compared to sim i lar studies of type-II QSOs 
dHumphrev et al.1 l20Tot iGreene et al.l I20TTI) . the NLR of 
the two type-I QSOs studied here appears to be less extended. 
Due to the small sample size, and the intrinsic difficulty of 
inferring the true "size" of emission in long-slit observations, 
we are unable to draw a conclusion on whether type-I QSOs 
actually have less extended NLRs. Recently, Husemann et aP 
(2013) carried out an integral-field spectroscopy of a large 
sample of nearby type-I QSOs revealing a typical size 
of 10 kpc, similar to that of type-II QSOs studied by 
iGreene et al.l (1201 Ih . They also caution that the difference 
between the size determined by long-slit spectroscopy and 
the size based on a two-dimensional spectroscopy is on 
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average a factor of two, and can be as high as a factor of six. 
Ev en so, the mean s i ze of N LRs in type-I QSOs determined 
by iHusemann et alj ( 120 13h seems to be smaller than that 
of type-IIs (28 k pc) based on a comparable integral field 
spectroscopy by iLiu et alJ (120131) . although these type-II 
QSOs are at a slightly higher redshift. 

The lack of galacti c scale outflows in typ e-I QSOs has also 
been pointed out bv iHusemann et all (120 13h . Out of 31 type-I 
QSOs in their study, only three show kpc-scale outflow, and 
the spectra of extended emission generally show narrower 
[O III] line profile than the central QSO spectrum. In_the 
framework of the unified model of AGNs ( Antonuccill 1 993h . 
while the outflow in the projected plane may look less ex- 
tended in type-Is than type-IIs, the projected velocity is ex- 
pected to be larger as we look directly into the ionization cone. 
However, it is not yet clear whether type-Is have larger pro- 
jected outflow velocities than type-IIs (see Hus emann et alJ 
[2013; Liu et al. 2013). At the same time, although the [O III] 
emission is expected to be more elongated or biconical for 
type-II QSOs, the morphology of [O III] emission seems to 
be more elongated in type-I QSOs ( IHusemann et al.ll2013h . It 
is possible that in a merger-driven evolutionary scenario (e.g., 
[Sanders et al. 1988), AGN feedback in the ionized gas occurs 
mainly in obscured (type-II) phase. The sample probed here 
may still not be quite luminous enough to drive galactic wide 
energetic outflows that expel gas from the host galaxy and 
quench star formation as those found at high-redshift (e.g., 
iCano-Diaz et alJ[2012h . Solid evidence for large-scale out- 
flows related to AGN feedback still remains elusive. 
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